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Supergiant Fast X-ray Transients (SFXTs) are a new class of High Mass X-ray Binaries, discov- 
ered by the INTEGRAL satellite, which display flares lasting from minutes to hours, with peak 
luminosity of 10^^ — 10^^ erg s^'. Outside the bright outbursts, they show a frequent long-term 
flaring activity reaching an X-ray luminosity level of 10^^ — 10^^ erg s^^ as recently observed 
with the Swift satellite. Since a few persistent High Mass X-ray Binaries (HMXBs) with super- 
giant donors show flares with properties similar to those observed in SFXTs, it has been suggested 
that the flaring activity in both classes could be produced by the same mechanism, probably the 
accretion of clumps composing the supergiant wind. We have developed a new clumpy wind 
model for OB supergiants with both a spherical and a non spherical symmetry for the outflow. We 
have investigated the effects of the accretion of a clumpy wind onto a neutron star in both classes 
of persistent and transient HMXBs. 
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1. Introduction 

In the last seven years, the hard X-ray INTEGRAL observatory discovered many new hard X- 
ray sources [jl]]. In particular, almost 30% of the new discovered sources are HMXBs, which were 
not detected in earUer observations. Among these, INTEGRAL discovered two classes of HMXBs 
with supergiant companions: the first class is composed of intrinsically highly absorbed hard X- 
ray sources (e.g. IGR J16318-4848) [|2|]. The members of the second class, called Supergiant 
Fast X-ray Transients (SFXTs; |Q, ^), exhibit outbursts with duration of a few days composed by 
many flares lasting from minutes to a few hours as discovered by [js), |^ with the Swift monitoring 
of 4 SFXTs (IGR J16479-4514, XTE J1739-302, IGR J17544-2619 and AX J1841.0-0536). 
The behaviour of SFXTs is characterized by a high dynamic range, spanning 3 to 5 orders of 
magnitude, from a quiescent state at 10^^ — 10^^ erg s^^ up to the peak luminosity during outbursts 
of 10^^ — 10^' erg s^^ Swift also discovered that SFXTs display a fainter flaring activity with 
luminosities of 10^^ — lO^'* erg s^^ 

Many different mechanisms have been suggested to explain the SFXT behaviour: ||] pro- 
posed that the high dynamic range shown by SFXTs is due to transitions across the neutron star 
centrifugal barrier produced by a change in the donor wind density. In particular, [§] proposed 
that what distinguishes SFXTs from persistent HMXBs with supergiant companions is that SFXTs 
host magnetars with large spin period (~ 10^ s). Another possibility involves the presence of an 
equatorial wind component denser than the polar wind, and inclined with respect to the orbital 
plane of the compact object. In this framework, the outburst is produced when the compact ob- 
ject crosses the equatorial wind component and, consequently, accretes more matter |^ [l^. This 
mechanism has been successfully applied to the SFXT IGR J 11215-5952, which shows periodic 



outbursts {Porh ~ 165 days, [^). [ pT| ] proposed that the flaring activity in SFXTs is due to the 
sudden accretion of dense blobs of matter composing the supergiant wind. In the framework of 
the clumpy wind model proposed by [|T^], [ p3[ ] suggested that different orbital separations could 
play a role in the different behaviour of SFXTs and persistent HMXBs. Persistent HMXBs have a 
small orbital period, with a distance supergiant-compact object < 2 stellar radii, while in SFXTs 
the compact object orbits the companion at larger distances. 



2. A new clumpy stellar wind model 

Recently we developed a new clumpy stellar wind model for OB supergiants in HMXBs []T4|]. 
Assuming that OB supergiants are surrounded by a clumpy and spherically symmetric wind, we 
assumed for the first time that clumps follow a power law mass distribution 

where M^i is the mass of the clump, and [Ma - Mb] is the mass range. The rate of clumps produced 
by the supergiant is related to the total mass loss rate Mtot by means of the k parameter, and we 
defined / = Mci/Mtot as the fraction of mass lost in clumps, where M^i is the component of mass 
loss rate due to the clumps. We assumed spherical clumps, with radii R^, then we also introduced 
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a power law distribution of radii, Rd '. 

A^Md °^ clumps s^^ (2.2) 
Spectroscopic observations of O stars suggest that clumps have, on average, the same velocity law 



of a smooth stellar wind [|15[]. We can then assume for the clump velocity profile Vci(r): 

vci(r) = v„o f 1-0.9983— y (2.3) 



where Voo is the terminal wind velocity, 0.9983 ensures that v(/?ob) 0, Rqb is the radius of the 
supergiant and j8 is a constant [|l6|]. 

From the balance pressure equation and the continuity equation [|l^, |l^ we find the law de- 
scribing how the clump size increases with the distance from the supergiant star: 



Rclir) =RciiRs) 



2 / M 1/3 

Rh(R.) 



(2.4) 



where R^ is the sonic radius, where the clumps start outflowing from the star JI^]. 

For any given mass of the clump, we derived the lower-limit for the clump radius, starting 
from the assumption that, in order to be accreted by the compact object, the clump must escape 
from the supergiant, i.e. the radiative force due to the scattering of the ions of the clump with the 
UV photons must dominate over the gravity of the supergiant. We also derived the upper-limit 
for the clump radius, starting from the definition of the clump as a density enhancement in the 
smooth stellar wind: clumps with radii larger than the upper-limit would be less dense than the 
smooth stellar wind (inter-clump medium), in contrast with the clump definition [|l4|]. The upper 
and lower-limits for the clump radius are drawn in Fig. |l|. 

To calculate the X-ray luminosity produced by the accretion of the inhomogeneous wind, we 
modified the Bondi-Hoyle accretion model. Assuming different orbital configurations and clumpy 
wind properties, we found that the observative characteristics of the flares (luminosity, duration, 
number of flares produced), do not depend only on the orbital parameters, but potentially are also 
significantly affected by the properties of the clumps. This model has been successfully applied 
to four HMXBs: Vela X-1, 4U 1700-377, IGR J18483-0311, and IGR Jl 1215-5952 Q, H^. For 
IGR J11215-5952 we were able to reproduce the lightcurve observed with Swift, with the intro- 
duction of a clumpy equatorial wind component around the supergiant. This result is in agreement 
with the accretion mechanism proposed by [P] for this source. 



3. Comparison with the HMXB 4U 1700-377 



4U 1700—377 [ ]21| ] is a bright eclipsing X-ray binary (Porb = 3.412 d) composed by a compact 
object (a neutron star or a black hole), and the 06.5 Iaf+ star HD 153919, located at a distance of 
1.9 kpc [p2|]. This source is characterized by a strong flaring activity with variations as large as a 



factor of 10 — 100 on short time scales (from minutes to hours) [23]. 

We analyzed the IBIS/ISGRI pubhc data archive from 2003 March 12 to 2003 April 22, and 
from 2004 February 2 to 2004 March 1, for a net exposure time of ~ 5.2 days (excluding the 
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Figure 1: Upper (dashed line) and lower-limit (solid line) for the clump radius at the sonic radius These 
functions have been obtained assuming the following parameters: Mqb — 30 Mq, /?ob = 23.8 R©, Mtot — 
Mq yr-i, v„ = 1700 km s^^ j3 = 1, Md/Mwind = 0.7. 



eclipse phase). The data reduction was carried out using OSA 7.0, and from the extracted light 
curve (15 — 60 keV) we found a total of 123 flares. For each flare we extracted the spectrum in the 
range 22 — 100 keV. All the spectra are well fitted by a thermal Comptonization model (COMPTT 
in XSPEC). For each flare, we derived the 1 — 200 keV luminosity, which is always greater than 
5.8 X 10^^ erg s^^ We then applied our clumpy wind model to the INTEGRAL observations of 
4U 1700—377. We first compared the observed distributions of the flare luminosities and durations 
with those computed adopting the system parameters found by [ 24 1 : the supergiant has a luminosity 
log(L/L0) = 5.82 lb 0.07, an effective temperature Teff « 35000 K, radius /?ob ~ 21.9 R©, mass 
^OB ~ 58 M©; the mass of the compact object is Mx = 2.44 M©. As shown in Figure ^, the 
flare properties are well reproduced with our clumpy wind model for Mtot = 1.3 x 10^^ M0 yr^^. 
Ma = 5 X 10^^ g and Mb = 2 X lO'"^ g,C, = \.l,y= -6.5 and / = 0.75. We found that the numbers 
of observed (123) and calculated flares (116) are in good agreement. 
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Figure 2: Comparison of the 4U 1700—377 flare luminosities and durations distributions, as observed with 
IBIS/ISGRI (solid line), with that calculated (dashed line), assuming the following binary system parameters: 
Mob =58 M©, Rq^ = 21.9 R©, Mns = 2.44 M©, 7?ns = 10 km. The parameters for the supergiant wind 
are: Mtot = 1-3 x 10"*^ Me, yr^^ Voo = 1700 km s"', j3 = 1.3, vq = 10 km s^^ Ma = 5 x lO'*^ g and 
Mb = 2 X 10''' g, C = 1 -2, 7 = -6.5 and / = 0.75. 
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